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Abstract 

A systematic first-principles study within density functional theory on the geometrical struc- 
tures and electronic properties of unconventional fullerene C64 and its derivatives C64X4 (X = 
H, F, CI, Br) has been performed. By searching through all 3465 isomers of C64, the ground state 
of C64 is found to be spherical shape with D2 symmetry, which differs from the parent cage of the 
recently synthesized C64H4 that is pear-shaped with Cat, symmetry. We found that the addition 
of the halogen atoms like F, CI, Br to the pentagon-pentagon fusion vertex of C64 cage could en- 
hance the stability, forming the unconventional fullerenes C64X4. The Mulliken charge populations, 
LUMO-HOMO gap energies and density of states are calculated, showing that different halogen 
atoms added to C64 will cause remarkably different charge populations of the C64X4 molecule; the 
chemical deriving could enlarge the energy gaps and affect the electronic structures distinctly. It 
is unveiled that C64F4 is even more stable than C64H4, as the C-X bond energy of the former 
is higher than that of the latter. The computed spectra of C64H4 molecules agree well with the 
experimental data; the IR, Raman, NMR spectra of C64X4 (X = F, CI, Br) are also calculated to 
stimulate further experimental investigations. Finally, it is uncovered by total energy calculations 
that C64X4 could form a stable hexagonal monolayer. 
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Fullerenes are closed-cage carbon molecules of 12 pentagons and several hexagons. Since 
discovery of Buckminster fuUerene Cgo H in 1985, intensive theoretical and experimental 
works have been done to study the structures and properties of fullerenes. The conventional 
fullerenes such as Cm, C70, Cjo+2n,n = 1,2,... are satisfying the empirical isolated pen- 
tagon rule (IPR) y, 0] , which states that the most stable fullerenes are those in which every 
pentagon is surrounded by five hexagons, ensuring the minimum curvature of the fullerene 
cage (i.e., the minimum strain) j4|. Fullerenes with other number of carbons are impossible 
to have all the pentagons isolated in the carbon cage, thus they are called unconventional 
fullerenes or non-IPR fullerenes, which are believed to be labile and difficult to synthesize. 

n 

In fact, most of the several stable fullerene isomers [5|] which had been isolated are satisfying 
IPR. For applications, fullerenes and their derivatives are potential new agents and mate- 
rials for molecular electronics, nanoprobes, superconductors, and nonlinear optics. Thus if 
IPR could be efficiently violated in some ways, the bank of stable fullerenes for practical 
applications will be much enriched. 

In the recent years, it has been uncovered that some non^IPR fullerenes may be stabilized 
through metallic endohedral (Sc2@C66[-0L A^Sc3_:i.N@C68[7], La2@C72 8], etc.). The other 
effort and progress have also been made in chemical deriving. In 2004, Xie et a/. |9!] had 
successfully synthesized a new exohedral chemical deriving non-IPR Ds^ fullerene [50] (say, 
C50CI10, with 10 chlorine atoms added to the pentagon-pentagon vertex fusions) in milligram 
by a graphite arc-discharge process modified by introducing a small amount of carbon tetra- 
chloride (CCI4) in the helium atmosphere. Later, another new chemical deriving non-IPR 
fullerene [64] (namely, C64H4, with 4 hydrogen atoms added to the vertex of a triplet directly 
fused pentagons) has been successfully prepared by Wang et al. 1^ in milligram with similar 
method. As CCI4 was added to provide chlorine source in the case of C50CI10, methane (CH4) 
was introduced to provide hydrogen for C64H4. Both C50CI10 and C64H4 are stabilized by the 
chemical modification on the pentagon-pentagon fusion vertex. This clue hints a promising 
general method to stabilize non-IPR fullerenes. While an efficient and controllable method 
to prepare the metallic endohedral fullerenes with macroscopic quantity is not yet matured, 
the chemical deriving might open a useful way to obtain stable non-IPR fullerenes. 

Moreover, many questions are still unclear: What is the ground state of C64? How about 
the electronic properties? Since the adding atoms in C50CI10 are CI, what will happen if 
the H atoms in C64H4 are replaced by other atoms, for instance, F, CI, Br? Here, by first- 
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principles density functional calculations, we shall address these issues, and report intensive 
studies on the geometrical structures, electronic properties and spectroscopies of and 
its derivatives C64X4 {X = H, F, CI, Br). As it is uncovered that C50CI10 molecules could 
form a hexagonal monolayer we also investigate the possibility that C64X4 could form 
a hexagonal monolayer. Our present study would deepen understanding of the properties of 
the recently isolated fulleride C64H4, and could also stimulate further experimental efforts 
on C64X4 (X = F,Cl,Br). 

We ha ve g enerated all of the 3465 topologically distinct isomers of C64 cage by spiral 
algorithmic], and sorted them by the corresponding pentagon-pentagon adjacency (PPA) 
numbers. (Definition: one distinct C-C bond shared by two adjacent pentagons in fullerenes 
is counted as one PPA.) Then, we have calculated the 6^4 isomers with small PPA numbers 
(AppA=2,3,4) by the semiempirical MNDO method[13]. (see supporting materials.) In 
Gaussian-like terminology, that would be called MNDO//MNDO calculations, in which 
B / / A represents B level energy calculation based on the geometry optimized on A level. It 
had been proved that the semiempirical MNDO method could reproduce the ab initio or 

nn 

density functional relative energies of fuUerene isomers with useful high accuracy |14J. I15|. 

Among them, nine isomers of C64 with the lowest energies (Fig. 1) were chosen, and 
their derivatives were g e,..ted. „s„, t.. density f„,.tio„a, theo.. (DFT)H, we have 
calculated geometries and relative energies of above nine isomers and their derivatives. In 
particular, we have invoked the B3LYP functional, which combines the Becke three pa- 



rameter non-loca 
of Lee- Yang- Parr 



hybrid exchange potential[l7t| and the nonlocal correlation functional 



181] . The basis set 6-31G* was used in the B3LYP calculations. The 



B3LYP/6-31G* calculations also give the LUMO-HOMO gaps. The harmonic vibrational 
frequencies and nuclear magnetic shielding tensors were calculated at the B3LYP/6-31G* 
and GIAO-B3LYP/6-31G* level, respectively. 

All of the calculations are completed ,,Ga„...03H except that the density of states 
(DOS) and monolayer total energies of C64X4 are calculated by means of the ABINIT 
package I20I. This package is coded within the DFT framework based on pseudopoten- 
tials and plane waves. Particularly, in our calculations, TrouUier- Mart ins norm conserv- 
ing pseudopotentials[21] and the Teter parametrization 22|] of the Ceperley-Alder exchange- 
correlation potential were used. 

Fig. 1 shows shapes and symmetries of the nine lowest-energy isomers of C64 obtained 
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3465:1 91 1(C3v) 3465j3402(Ci) 3465:3423(C, 




3465:3425(Ci) 3465:3451 (D2) 3465:3452(Cs) 




FIG. 1: The 9 lowest-energy isomers of C64 obtained at MNDO. 



at MNDO level. They have been labeled as 3465 : n , where n denotes the isomer number 
in spiral nomenclature [l^. We reoptimized the geometries and recalculated their relative 
energies by B3LYP/6-31G* method, where the LUMO-HOMO gap energies are obtained. 
The details are listed in Table [B 



As seen in Table [U the isomer 3451 (numbering in spiral nomenclature 12]) is the ground 
state of C64, followed by the isomer 3452 and 3457. The isomer 3451 is nearly spherical 
shaped, with D2 symmetry, while the isomer 3452 and 3457 bear Cg and C2 symmetry, re- 
spectively. These three lowest-energy isomers have two pairs of adjacent pentagons (counted 
as 2-PPA), and the rest isomers with higher energies all have three couples of adjacent pen- 
tagons (counted as 3-PPA) except that the isomer 1911, has a triplet of fused pentagons 
(also counted as 3-PPA). This pear-shaped Ca^ symmetry isomer 1911, which is just the 
carbon cage parent of the synthesized C64H4 by Wang et a/. [lOj, however, is not apparently 
favored in energy when it is compared with other eight isomers. The ground state isomer 
3451 also has the largest LUMO-HOMO gap energy, followed by the isomer 3457, whereas 
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TABLE I: MNDO//MNDO, B3LYP/6-31G*//B3LYP/6-31G* relative energies (in eV), LUMO- 
HOMO gap energies (Egap, in eV) of the nine lowest -energy Cg4 isomers 

Relative Energy 



Isomer'^ Symm NppA^ Npppv' MNDO'^ B3LYP'^ E, 
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(1911) 
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1.0979 


[4] 


1.2403 


[6] 


1.0879 
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1 4144 
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1 4941^ 
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(3423) 


Cs 
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1.1597 


[6] 


1.2156 


[5] 


1.3845 


4 


(3425) 


Ci 


3 


6 


1.4319 


[9] 


1.2995 


[7] 


1.2408 


5 


(3451) 


D2 


2 


4 


0.0000 


[1] 


0.0000 


[1] 


2.2327 


6 


(3452) 


Cs 


2 


4 


0.2600 


[2] 


0.2987 


[2] 


2.0482 
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(3454) 


C2 


3 


6 


1.3619 


[7] 


1.4330 


[8] 


1.6311 


8 


(3455) 


Ci 


3 


6 


1.1178 


[5] 


1.1116 


[4] 


1.9369 


9 


(3457) 


Cs 


2 


4 


0.4807 


[3] 


0.5563 


[3] 


1.8547 



"Labeled as in(n) form, where (n) was the numbering in spiral nomenclature 

^Number of pentagon-pentagon adjacency (PPA). 

'^Number of pentagon-pentagon fusion vertex (PPFV). 

''[n] indicates the stability order on MNDO level. 

■^[n] indicates the stability order on B3LYP/6-31G* level. 

that of the isomer 1911 is the smallest in the nine. 

It may be worth to note that the stability order of the nine isomers obtained by B3LYP/6- 
31G* calculations is similar, but not exactly the same to that by MNDO. Both methods give 
the same order for the most stable structures (the isomer 3451, 3452, 3457). 

The pentagon-pentagon fusion vertex (PPFV) are usually considered to be the most 
active sites of the fullerenej23|, and the syntheses of stable C50CI10 and C64H4 also clearly 
evidence that the chemical modification on the PPFV could effectively release the large steric 
strain of non-IPR fuUerenes. Here we attach H and CI atoms to the PPFV of above nine C64 
isomers to construct C64 fuUerene derivatives. As listed in TablelU isomers 1911, 3451, 3452, 
3457 have 4 PPFV, while others have 6 PPFV. Thus, the corresponding fullerene derivatives 
are C64X4 and 064X5 {X = H or CI), respectively. Since the number of the added X atom 
is not fixed, the relative energy is not suitable for evaluating the relative stability. In order 
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TABLE II: B3LYP/6-31GV/B3LYP/6-31G* bond energies AE (in eV), LUMO-HOMO gap ener- 
gies (Egap , in eV) of the derivatives of nine lowest-energy C64 isomers 
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2 ('3402) 


6 


3 1113f4l 


0.9611 


1 8532 f4l 


8963 


3(3423) 


6 


3 0093 f7l 
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2.9186[8] 


0.6833 


1.7008 [8] 


0.7959 


5(3451) 
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2.8985 [9] 


1.3557 


1.6392 [9] 


1.2398 


6(3452) 
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3.0153[6] 


1.3968 


1.7628[6] 


1.3685 


7(3454) 


6 


3.1407[3] 


1.7949 


1.8827[3] 


1.7374 


8(3455) 


6 


3.0754[5] 


1.2651 


1.8214[5] 


1.2441 


9(3457) 


4 


3.1920[2] 


1.9108 


1.9412[2] 


1.8623 



"the number of added H or CI atoms, i. e., the number of pentagon-pentagon fusion vertex (PPFV) for 
each isomer according to the method we generated the derivatives. 

to discuss the relative stability of Cq^Xu, we define C-X bond energy AE as follows |2^: 

AE = ^{Ec,, + nE^~Ec,,y,J 

The higher the C-X bond energy AE, the greater the stability of C64X„. 

Table [TTl hsts the B3LYP/6-31G7/B3LYP/6-31G* C-X bond energy AE and relative 
stability, LUMO-HOMO gap energies of the hydro- and chloro- derivatives of nine lowest- 
energy C64 isomers such as C64X„ {X = H or CI, n = 4 or 6). We still use the spiral 
nomenclature to label the generated fullerene derivatives. 

The stability order is changed after H atoms are attached to the C64 isomers. As expected, 
the C64H4 1911, with a remarkably largest C-H bond energy, turns out to be the most stable 
in all C64H„ derivatives. The following ones are C64H4 3457 and C64H6 3454, both having C2 
symmetry. Derivative 3451, which has been evolved from the ground state of C64 isomers, 
becomes the most instable among the nine derivatives. 

Corresponding to its remarkable larger C-X bond energy AE, the C64H4 1911 has the 
largest LUMO-HOMO gap energy, 2.6697 eV, while that of its parent C64 1911 is only 
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1.0879 eV. The LUMO-HOMO gap energies of others are also shifted, implying that the 
chemical deriving on PPFV could affect the electronic structures of non-IPR fuUerenes. 

Obviously, the chemical deriving on PPFV could indeed enhance the stability of the non- 
IPR fullerenes in the case from 1911 to C64H4 1911. However, the effects are different 
for different isomers. It seems that the vertex of the triplet fused pentagons are rather more 
favorable for hydrogen addition than other types of vertex, which is consistent with the 



suggestion by Liu et al. 



25|. 



The effects of CI addition are similar to that of the H addition. After CI atoms attached 
to the C64 isomers, the stability order is changed, and the new order is entirely the same as 
that of the H addition. C64CI4 1911 also becomes the most stable with a remarkable larger 
C-Cl bond energy. The LUMO-HOMO gap energies are changed with the same manner 
too. C64CI4 1911 also has the highest LUMO-HOMO gap energy, 2.6923 eV, in all C64C1„ 
derivatives. These give us a clue that haloid atoms may be appropriate for the chemical 
deriving on the PPFV of non-IPR fullerenes to stabilize the structure. 

Now we attach H, F, CI, Br atoms to the PPFV of the isomer C64 1911 to construct 
C64X4 {X = H, F,Cl, Br) fullerene derivatives and, to study their structural and electronic 
properties. Fig. 2 displays the geometrical structures of C64X4 (X = H, F,Cl, Br). These 
molecules like pears in shape with four short hairs on its tapering top or odd pears with four 
stalks. The tapering top of C64 cage is composed of a triplet of directly fused pentagons. 
Its strain is obviously high due to its high surface curvature, which is a key factor for its 
instability. The length of C-X bond along the C3 axis in C64X4 {X = H, F,Cl, Br) is about 
1.09, 1.36, 1.77, 1.93 A, respectively, while the length of other C-X (X = H,F,Cl,Br) 
bonds around the C3 axis is about 1.10, 1.38, 1.82, 1.99 A, respectively. The bond lengths 
of C-X bonds along the C3 axis are nearly the same as the experimental value of C-X 
{X = H, CI, Br) bond length in methane (CH4), CF4, chloroform (CCI4) or CBr4, while the 
other bonds around the C3 axis are somewhat longer. For a comparison, we have tried to 
put H, F, CI or Br atom along the C3 axis into the inner of C64 cage, and found that the 
calculated energies are higher than those when these atoms are added to the outer of the 
cage, suggesting that the exohedral fullerenes C64X4 {X = H, F, CI, Br) are more stable. 

Fig. 2 also show the MuUiken atomic charge populations of C64X4 (X = H, F,Cl, Br). 
The PPFV C atoms (i.e. those bonded with X atoms) in C64H4 are colored in red, indicating 
that their charges are rather negative. The other atoms at the tapering top of 6*64 cage 
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FIG. 2: Schematic structures and Mulliken atomic charge populations of Csv C64X4 molecules: (a) 
C64H4, (b) C64F4, (c) C64CI4, and (d) C64Br4. The capital letters indicate the exohedral atoms, 
and the rest atoms are all C atoms. The color of the atoms show the Mulliken atomic charge 
amount in electrons. 

are marked in dark green, showing that they carry little positive charges. The atoms at 
bottom are electroneutral as they are in black color. The charge populations of C64X4 
(X = CI, Br) are mainly similar to C64H4. The PPFV C atoms in 054X4 (X = CI, Br) are 
also red, negative charged; while the tapering top of Cq^ cage are dark green, little positive 
charged. Interestingly, the colors of X atoms (namely light green, dark green, and dark 
red for X = H,Cl,Br, respectively), are quite distinct, revealing that the charge transfers 
of H, CI, Br atoms, when fused to the 6*64 cage, are different. The charges of H atoms 
are rather positive, while those of CI atoms are less positive, and those of Br atoms are 
even little negative. With regard to electron transfers, the PPFV C atoms act as electron 
acceptors, while H and CI atoms act as electron donors, but Br could not. Surprisingly, the 
addition of F atoms behaves rather differently. Contrary to C64X4 (X = H,Cl, Br), the 
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PPFV C atoms are positive charged (bright green color) and F atoms are negative charged 
(hght red color), showing a totally inverse electron transfers. All the above results imply 
that a proper chemical modification on PPFV could change the electronic structures of the 
non-IPR carbon cage, but adding different atoms would bring different consequences. 

For C64H4 and C64CI4, we have calculated the C-X bond energy AE. In the same method, 
the obtained C-X bond energy AE for C64X4 {X = F, Br) are 4.2810 eV, 2.4200 eV, respec- 
tively. As AE of C64H4 is 2.6697 eV, it seems that C64F4 is even more stable than C64H4. 

We have also calculated the LUMO-HOMO gap of C64 and C64X4 {X = H, F, CI, Br) by 
B3LYP/6-31G*. The gap of C64 is 1.0879 eV, while those of C64X4 are 2.6697 eV, 2.7075 
eV, 2.6923 eV, 2.6762 eV for X = H, F, CI, Br, respectively, which is comparable to the gap 
energy (2.758eV) of Cqo calculated within the same computational method. A large LUMO- 
HOMO gap energy has been taken as an indication of a large stability of fullerenes 26|, H]- 



The remarkable enlargement of the gap energy from C64 to C64X4 implies that the chemical 
deriving could indeed enhance the stability of non-IPR fullerenes. 

The DOS provide a convenient overall view of the electronic structure of clusters and 
solids. Fig. 3 exhibits the DOS of C64 and C64X4 {X = H, F,Cl, Br), which are obtained 
by Fermi-Dirac smearing of molecular energy levels. We can see that from C64 to C64X4 
(X = H, F, CI, Br), LUMO-HOMO gaps are enlarged. 

To provide a verifying basis for the experimental identification of C64X4 (X = 
H, F,Cl, Br), we have calculated the IR, Raman spectra by B3LYP/6-31G*, and NMR 
spectra by GIAO-B3LYP/6-31G*, with the optimized geometry in the B3LYP/6-31G* level, 
as shown in Figs. 4 and 5. 

As can be seen in Figs. 4(a) and 4(b), the computed IR spectrum of C64H4 agrees well 
with the experimental resultsflOl]; both IR and Raman spectra of C64H4 spread up to more 
than 3000 cm~^ and have two major regions. The region of two sharp peaks near 3000 cm~^ 
corresponds to the C-H stretching modes which are rather active in both IR and Raman 
spectra, and the wide region in the range less than 1700 cm^^ corresponds to the C-C 
stretching, C-C-C bending and C-C-H bending modes in which active modes are different 
for IR and Raman spectra. 

The IR and Raman spectra of C64X4 {X = F,Cl, Br) are rather different, as depicted 
in Figs. 4(c) to 4(h). Unlike C64H4, they are only extended to circa 1700 cm~^ and have 
no distinct separated regions. From a view of vibrational modes, the spectra of C64X4 
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FIG. 3: Density of states (DOS) and the integrated DOS(doted line) of electrons as a function of 
energy for C64 and C64X4 {X = H,F,Cl,Br). The energy zero-point is taken at the Fermi level, 
that is indicated by the vertical dash-dot line. 

{X — F, CI, Br) lack a clear region of C-X stretching modes like C64H4. It may be because 
H atom is much lighter than C atom, the C-H stretching modes could be separated from 
other modes caused by C atoms, and the corresponding peaks are located at a very high 
frequency region. By contrast, F, CI and Br atoms are heavier than C atom, the motion of 
F, CI and Br atoms must be coupled with C atoms strongly, resulting in that the IR and 
Raman spectra of C64CI4 could not be parted into two regions. This may also cause other 
spectral differences between C64H4 and C64CI4. 

For C64X4 {X — F,Cl,Br), the IR spectra profiles exhibit some similarity, but the 
main peaks are located at different frequencies. It seems that, from CgiFi to C64Br4, the 
main peaks of IR spectra are red shifted, and intensities are reduced. Raman spectra show 
more similarity: not only are the profiles very similar, but also the main peaks are all 
located between 1400 cm~^ and 1600 cm~^. By a careful check, the Raman spectrum of 
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C64H4 below 1600 cm~^ in Fig. 4(b) also presents the above similarity to that of C64X4 
{X — F, CI, Br). In addition, the Raman spectral activity grows from C64H4 to C64Br4. We 
can conclude that the IR spectra are more sensitive to different added X atoms than Raman 
spectra: the heavier added X atoms will red-shift the main active IR modes and reduce their 
intensities, but will enhance the activity of Raman modes and not affect the active Raman 
modes frequencies. 

Fig. 5 gives the calculated NMR spectra of C64X4 {X — H, F,Cl, Br). The computed 
NMR spectrum of C64H4 agrees well with the experimental results as depicted in Fig. 4(a). 
The Csv C64X4 have 14 unique types of C atoms [60 sp^ C atoms (8 x 6; 4 x 3); 4 sp^ C atoms 
(1 X 3; 1 X 1)] and 2 unique types of X atoms [total 4 X atoms (1 x 3; 1 x 1)], where a x b 
represents a unique types, and in one type, there arc b symmetrical equivalent atoms. All of 
the calculated C^^ NMR spectra of C atoms in C64X4 {X — H, F, CI, Br) are easily divided 
into two regions. One may see that there are 12 peaks in the left region, 8 of which have 
the same intensity, and the other 4 have the half of the intensity, which could be ascribed 
to the sp^(8 x 6; 4 x 3) C atoms. The right regions of C64X4 {X = if, F, CI, Br) only have 
two peaks, whose relative intensities are 3:1. It is interesting to note that the two peaks in 
C64H4 are very close, where the higher peak is in the left, whereas the two peaks in C64F4 
are separated from each other, where the higher peak is in the right. We can ascribe this 
feature to the sp^(l x 3; 1 x 1) C atoms which are bonded to the X atoms (i.e. PPFV). The 
lower peak could be ascribed to the single C atom on the C3 axis, and another peak could 
be ascribed to the equivalent triple C atoms out of the C3 axis. The NMR spectra of C64CI4 
and C64Br4 are rather similar to those of C64F4, with small quantitative differences. The 
calculated NMR spectrum of H atoms in C64H4 is also included, where the two peaks 
are observed in an intensity ratio of about 3:1. Obviously, the lower peak can be ascribed to 
the single H atom on the C3 axis, whereas another peak can be ascribed to the equivalent 
triple H atoms out of the C3 axis. 

Let us now address another interesting issue if C64X4 [X — H, F, CI, Br) molecules can 
condense to form a monolayer structure. We have considered the hexagonal lattice, which 
has been simulated by a supercell with layer-layer distance of 30 A. The C64X4 molecules are 
put on the lattice site with C3 axis perpendicular to the monolayer plane. Then, the total 
energy is calculated accordingly with respect to different lattice parameters (every 0.1 A) for 
the presumed structure. For C64H4 hexagonal monolayer, as indicated in Fig. 6, a clearly 
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single minimal total energy with the optimized lattice constant ao = 9.89 A is obtained, 
revealing that this structure is stable and easily to form. The data with 0.01 A precision 
are obtained by a sphne interpolation from the calculated data with 0.1 ^4 precision. The 
cohesive energy per C64H4 molecule for this monolayer is found to be 0.48 eV. For C64X4 
{X = F,Cl,Br), the hexagonal monolayer structure is also found to be stable, and the 
optimized lattice constants are still Oq = 9.89 A. The cohesive energy per molecule are 0.67 
eV, 0.55 eV, 0.32 eV, respectively. It appears that C64F4 hexagonal monolayer phase is the 
most stable among them. 

To summarize, in terms of first-principles calculations based on the density functional 
theory, we have systematically reported the structures, electronic properties and spectro- 
scopies of the unconventional fuUerene C64 and its derivatives C64X4 {X = H, F,Cl, Br). 
We have found that the parent carbon cage of the recently synthesized C64H4 with C^y 
symmetry is not the ground state of C64, while the latter should be a spherical shaped, D2 
symmetrical isomer with spiral number 3451. Our computation results show that the com- 
bination of the four X (X = H,Cl,F, Br) atoms to the pentagon-pentagon fusion vertex of 
C64 cage could change the electronic structures and enhance the stability of the fullerenes, 
and C64X4 {X = F, CI, Br) may be also stable unconventional fuUerides like C64H4. Similar 
to C64H4, C64X4 {X = F, CI, Br) are found to be pear-shaped, with large LUMO-HOMO 
gap energies and large C-X bond energy. It appears that C64F4 is even more stable than 
C64H4 as the C-X bond energy of the former is higher than that of the latter. The computed 
spectroscopies of C64H4 molecule agree well with the experimental data; the IR, Raman, 
NMR spectra of C64X4 {X = F, CI, Br) are also presented to stimulate further experimental 
investigations. In addition, it is revealed that, by total energy calculations, C64X4 could 
form stable hexagonal monolayers and among them, C64F4 hexagonal monolayer is the most 
stable. 
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FIG. 4: IR and Raman spectra of Csv C64X4 (X=H, F, CI, Br), which are scaled by a factor of 
0.9613. 
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FIG. 5: NMR spectra of Csv C64X4 {X = H, F, CI, Br), scaled by referenced shifts, where the inset 
of (a) shows the chemical shifts of H atoms, others show that of C atoms. 
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FIG. 6: The total energy per C64H4 molecule in hexagonal monolayer, measured from the total 

energy of the isolated C64H4 molecule, as a function of the lattice constant, where the optimized 
lattice constant, ao =9.89 A. 
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